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EXECUTIVE SUMMARY 


This report presents the results of recent investigations into the aerodynamics of simulated 
runback ice accretion on airfoils. Aerodynamic tests were performed on a full-scale model using 
a high-fidelity, ice-casting simulation at near-flight Reynolds (Re) number. In addition, follow- 
on subscale tests were conducted with low-fidelity simulations on a quarter-scale model at low 
Reynolds number. The ice-casting simulation was attached to the leading edge of a 72-inch 
(1828.8-mm) chord National Advisory Committee for Aeronautics (NACA) 23012 airfoil model. 
Aerodynamic performance tests were conducted at the Office National d’Etudes et de 
Recherches Aerospatiales (French Aeronautics and Space Research Center) (ONERA) FI 
pressurized wind tunnel over a Reynolds number range of 4.7xl0 6 to 16.0xl0 6 and a Mach (M) 
number range of 0.10 to 0.28. For Re = 16.0 X 10 6 and M = 0.20, the simulated runback ice 
accretion on the airfoil decreased the maximum lift coefficient from 1.82 to 1.51 and decreased 
the stalling angle of attack from 18.1° to 15.0°. The pitching-moment slope was also increased 
and the drag coefficient was increased by more than a factor of two. In general, the performance 
effects were insensitive to Reynolds and Mach number changes over the range tested. 

The results from the full-scale tests were used to evaluate aerodynamic simulation methods for 
runback ice accretion. Aerodynamic tests were conducted on a quarter-scale NACA 23012 
model (18-inch (457.2-mm) chord) at Re = 1.8 X 10 6 and M = 0.18, using low-fidelity, 
geometrically scaled simulations of the full-scale casting. Simple, two-dimensional simulations 
of the upper- and lower-surface runback ridges provided the best representation of the full-scale, 
high Reynolds number iced-airfoil aerodynamics. Higher-fidelity simulations of the runback ice 
accretion that included geometrically scaled three-dimensional features resulted in larger 
performance degradations than the full-scale model. At this time, it is not clear why a simple, 
two-dimensional, geometrically scaled ridge simulation is required to best represent the full- 
scale, high Reynolds number, iced-airfoil aerodynamics. Reynolds number effects between the 
subscale model (1.8 X 10 6 ) and full-scale model (at 4.7 X 10 6 ) may be significant and more research 
is recommended. 

A new subclassification of spanwise ridge ice that distinguishes between short and tall ridges is 
proposed. This subclassification is based upon the flow field and resulting aerodynamic 
characteristics, regardless of the physical size of the ridge and the ice-accretion mechanism. Tall 
spanwise ridges have a profound effect on the airfoil flow field with a large (and often unsteady) 
separation bubble that grows rapidly with angle of attack and precipitates stall at a much lower 
lift coefficient and angle of attack than other ice accretions. In contrast, short spanwise ridges 
are characterized by a small, stable separation bubble formed in the immediate vicinity of the 
ridge. This small separation zone results in a limited effect on the airfoil pressure distribution 
relative to the clean configuration. Significant differences in surface pressure are only observed 
at the ridge location. The separation bubble does not significantly increase in size with angle of 
attack. Thus, the stall of the iced airfoil is generally governed by trading-edge separation 
moving forward with increasing angle of attack (trading-edge stall). In some cases, the airfoil 
may stall from the leading edge or from the ridge itself. More research is needed to determine 
the appropriate simulation methods for short ridges, particularly for Reynolds numbers less than 
1.8xl0 6 . 


vii 
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1. INTRODUCTION. 


Many aircraft certified for flight in icing conditions employ anti-icing systems on wing and 
empennage leading edges, as well as engine nacelles and other aerodynamic surfaces. In some 
cases, these ice-protection systems allow some fraction of the impinging water to flow 
downstream from the protected region and freeze, leaving a ridge-type ice accretion. 
Understanding the aerodynamic effects associated with runback icing is important to ensuring 
safety of flight operations and reducing development and certification costs. It is difficult and 
expensive to document the aerodynamic effect of this type of ice accretion in natural icing 
conditions. Therefore, simulations are often used to represent the ice accretion in a wind tunnel, 
sometimes on a subscale model at lower-than-flight Reynolds number. Subsequent questions 
naturally arise regarding the fidelity of the simulated ice accretion and Reynolds number effects. 

Many of these same issues surround iced-airfoil aerodynamics in general and have recently been 
addressed in the Airfoil Ice Accretion Aerodynamics Simulation research program sponsored by 
the National Aeronautics and Space Administration (NASA) and the Office National d’Etudes et 
de Recherches Aerospatiales (French Aeronautics and Space Research Center) (ONERA) [1]. 
This program quantified many of the uncertainties associated with various subscale simulation 
methods developed to capture the essential aerodynamic features of the full-scale iced airfoil. 
The authors considered four types of ice accretion classified by their unique aerodynamics: 
roughness, streamwise ice, horn ice, and spanwise ridge ice. The aerodynamic characteristics of 
the spanwise ridge ice accretion were developed from research related to large ice accretions 
downstream of protected surfaces typically formed in supercooled large drop (SLD) conditions. 
These may be referenced to as tall spanwise ridges, such that the flow field is characterized by a 
large separation bubble aft of the ridge. This leads to significant airfoil performance penalties, 
often surpassing that of horn ice. However, as shown in this report, runback ridges do not 
usually exhibit this type of flow field. Spanwise ridges resulting from anti-icing systems tend to 
be shorter in height above the airfoil surface, cover a larger surface length, and have a more 
streamlined geometry. They may be referenced to as “short” due to their reduced impact on the 
iced-airfoil aerodynamics. Regardless of the ice accretion mechanism, the aerodynamics of 
“tall” and “short” ridges are significantly different, thus leading to a subclassification of 
spanwise ridge ice. Bragg, et al. [2] allude to this case in their review of iced-airfoil 
aerodynamics. This comparison is covered in more detail in this report. 

Aerodynamic perfonnance data for runback ice accretion (or any high-fidelity simulation) on 
airfoils or wings are scarce in the public domain. Gray and von Glahn [3] documented the 
associated drag increases due to runback ice accretion in an icing wind tunnel, but were unable to 
conduct these and other perfonnance measurements over a significant angle of attack range. 
Calay, et al. [4] simulated runback ice accretion ridges using a step, a ramp, and a triangular 
shape, each with a normalized height, k/c = 0.0035 on a National Advisory Committee for 
Aeronautics (NACA) 0012 airfoil at Re = 1.25xl0 6 and M = 0.08. They found that the drag 
coefficient at 0°-angle of attack increased by up to approximately 0.011, while maximum lift 
coefficient decreased by up to 0.35 with the shapes at x/c = 0.05. Calay, et al. [4] also found that 
the same simulations, when located at x/c = 0.15, increased Ci, max of the airfoil. The greatest 
Ci,max increase was approximately 0.08 and a sta ii was delayed by 1°. The authors noted that the 
stall of the airfoil with the ice simulations began from the artificial ice shapes rather than from 
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the trailing edge, as was the case for the clean airfoil. The increase in Ci, max was attributed to the 
flow remaining attached at greater angles of attack than in the clean case. The mechanism by 
which this occurred was not discussed except to note that the shape added “extra turbulence.” 
Calay, et al. [4] concluded that small changes in the artificial ice shapes were able to produce 
large changes in performance effects, thus requiring accurate simulations to properly estimate 
aerodynamic effects. Papadakis and Gile-Laflin [5] also observed performance increases near 
maximum lift due to a backward facing ramp with k/c = 0.0041 at x/c = 0.15 and a spoiler with 
k/c = 0.0053 at x/c = 0.15. Their tests were conducted using a modified NACA 63a-213 airfoil at 
Re = 2. Ox 10 6 and M = 0. 17. The ramp increased Ci jmax by 0. 1 1 and delayed stall by 4° while the 
spoiler increased Ci, max by 0.01 and delayed stall by 1°. Tests with the ramp farther forward on 
the airfoil, at x/c = 0.025 reduced Ci, max by 0.23 and the stalling angle of attack by 2°. These 
reports indicate mixed results in perfonnance for simulated runback ridges. Relative to the 
clean, or uniced airfoil, increases in drag were always observed. However, these artificial ice 
shapes increased the maximum lift coefficient and stalling angle of attack in some cases. Such 
increases run contrary to the recent body of research associated with tall, SLD-type spanwise 
ridge ice accretion [6-8]. 

Whalen, et al. [9 and 10] investigated the aerodynamic perfonnance penalties of scaled, runback- 
type ice accretions on NACA 23012 and NACA 3415 airfoils at Re = 1.8xl0 6 and M = 0.18. 
The ice accretions that these simulations were based on were generated on a full-scale model at 
the NASA Glenn Icing Research Tunnel (IRT). The results indicated that runback ice accretion 
can significantly degrade the aerodynamic performance of airfoils, particularly for ridges located 
near the airfoil leading edge. It was also shown that, for the shortest accretions, these effects 
were more sensitive to the height of the simulation. 

Whalen, et al. [9 and 10] also found that the performance penalties associated with simple, two- 
dimensional simulations of runback ice accretion did not agree well with higher-fidelity 
simulations that included roughness and a more accurate ice-shape profile. Both simulations 
were geometrically scaled based upon the ratio of the chord length of the aerodynamic model to 
that of the icing model. Results similar to the aforementioned studies [4 and 5] were observed. 
A simple, square-cylinder simulation ( k/c = 0.0035) caused an increase in maximum lift and 
stalling angle of attack for the NACA 3415 airfoil and had a negligible effect on the maximum 
lift and stalling angle of attack of the NACA 23012. Further investigation revealed that this 
simulation was similar in height to the boundary-layer thickness. The source of this performance 
effect was linked to the shape’s ability to generate a mixing layer that entrained higher- 
momentum fluid into the boundary layer while generating a small, stable separation bubble. The 
effect was more pronounced in the case of the NACA 3415 because it exhibits a trailing-edge 
stall, unlike the NACA 23012, which stalls from the leading edge at the test Reynolds number. 
Trailing-edge stall can be effectively mitigated by increasing the momentum in the airfoil 
boundary layer. 

The effects observed for the runback-ridge simulations aroused suspicion as to the appropriate 
method for developing subscale simulations. Geometric scaling had been used in the past with 
good success [8 and 11]. Busch, et al. [12], in particular, provide quantitative assessment of ice 
accretion simulation methods on a quarter-scale model at low Reynolds numbers. The lift 
performance increases due to the two-dimensional simulations observed by Whalen, et al. [9 and 
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10] were generally not expected and did not agree with the results for higher-fidelity simulations. 
Because the simulations were similar in height to the local boundary-layer thickness, it was 
theorized that the boundary-layer thickness may be the appropriate length scale. Tests conducted 
by Whalen, et al. [9 and 10] using two-dimensional boundary-layer-scaled simulations, which 
were approximately twice as tall as the geometrically scaled simulations, showed significant 
penalties that were in better agreement with (geometrically scaled) higher-fidelity simulations of 
warm-hold runback ice accretion. However, without aerodynamic performance results at full- 
scale Reynolds number, it was difficult to draw a definitive conclusion concerning the scaling of 
these shapes. Lee, et al. [13] investigated geometry and Reynolds number scaling of a runback 
ice accretion on a business-jet wing. They found that geometric scaling did not reproduce the 
aerodynamics of the full-scale wing with the runback-type spanwise ridge. Empirical methods 
were used to develop aerodynamically equivalent shapes that represented the full-scale, iced- 
wing aerodynamics on the subscale wing. Reynolds number effects were also found to be 
significant for this ice shape, unlike the case of the large, leading-edge ice shape that was also 
tested. 

These recent results from subscale model tests at lower-than-flight Reynolds number clearly 
show the need for high-fidelity, high Reynolds number aerodynamic data. Therefore, a key 
objective of this work was to determine the effect of a high-fidelity runback ice accretion 
simulation on airfoil aerodynamic perfonnance at near-flight Reynolds number and quantify the 
Reynolds number effects. In addition, follow-on subscale tests were conducted with low-fidelity 
simulations at low Reynolds number to provide information about the accuracy of aerodynamic 
simulation methods for small runback-type ice accretion. The first objective was achieved by 
conducting aerodynamic performance tests at the ONERA FI pressurized wind-tunnel using a 
72-inch (1828.8-mm) chord NACA 23012 airfoil over a Reynolds number range of 4.7xl0 6 to 
16.0xl0 6 and a Mach number range of 0.10 to 0.28. The high-fidelity simulation of the runback 
ice accretion was developed from an anti-icing system test in the NASA IRT. The second 
objective was achieved by conducting aerodynamic performance tests in the University of 
Illinois low-speed wind tunnel using an 18-inch (457.2-mm) chord NACA 23012 airfoil at a 
Reynolds number of 1.8xl0 (1 and a Mach number of 0.18. Subscale simulations of the runback 
ice accretion, fabricated using simple-geometric shapes and roughness, were tested on the airfoil 
for comparison to the full-scale results. Finally, these results were used to further define the 
aerodynamic characteristics of spanwise ridge ice. A subclassification is proposed that 
distinguishes between tall and short spanwise ridges based upon their respective aerodynamic 
effects. 

2. EXPERIMENTAL METHODS . 

2. 1 THE ONERA FI WIND TUNNEL AND MODEL . 

The full-scale aerodynamic tests were performed at the ONERA FI pressurized wind tunnel 
facility, shown schematically in figure 1 [14]. The closed-return wind tunnel has a test section 
measuring 138 inches (3500 mm) high by 177 inches (4500 mm) wide by 433 inches 
(11000 mm) long. The maximum test section Mach number is 0.36 and the maximum stagnation 
pressure is 57 psia (3.85 bar). The unit Reynolds number can be varied up to a maximum of 
6.0xl0 6 /ft at Mach = 0.23. Total temperature is maintained via a heat exchanger (i.e., “water 
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cooler” in figure 1) located in the second diffuser downstream of the fan. The fan operates at 
constant speed while the test section Mach number is controlled by adjusting the pitch of the 
blades. The test section inlet flow is conditioned through a 7.18-to-l contraction containing a 
honeycomb flow straightener and three turbulence reduction screens. 


72 m 



Figure 1. Schematic of the ONERA FI Pressurized Wind Tunnel 

The 72-inch (1828.8-mm) chord NACA 23012 airfoil model was mounted vertically in the test 
section, as shown in figure 2. The model span was 137.48 inches (3492 mm) and was mounted 
in the floor force balance. Small gaps between the bottom of the model and the test section floor, 
as well as the top of the model and the test section ceiling, were maintained so as not to cause 
mechanical hysteresis in the force-balance measurements. The model had a main chordwise row 
of 72 pressure taps located at 43% span measured from the test section floor. In addition, there 
was a row of 20 taps oriented spanwise at xic = 0.70 on the upper surface. The model was 
designed and built with full-span removable, interchangeable leading-edge sections. The 
baseline leading edge had the clean NACA 23012 profile, while the alternate leading edge had a 
truncated nose geometry. The latter design facilitated mounting of the runback-ridge ice-casting 
simulation. Accommodations were also made in the pressure tapping to allow for rapid 
connection of pressure instrumentation in the cast ice shapes. Also shown in figure 2 is the wake 
rake located one chord length downstream of the model trailing edge. The wake rake had 100 
stagnation pressure probes spaced 0.79-inch (20-mm) apart and was located at a fixed spanwise 
station at 57% span above the test section floor. 
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Figure 2. Photograph of the Full-Scale NACA 23012 Airfoil Model Installed in ONERA FI 

Wind Tunnel Test Section 

Data acquisition runs were performed in angle of attack sweeps for increasing and then 
decreasing angle of attack at a constant sweep rate of 0.1°/sec. Data were also acquired at fixed 
angles of attack for selected angles over the range of the sweep, and repeat runs were perfonned. 
The data shown in this report are for increasing angle of attack sweeps and have been averaged 
to the nearest 0.5° in postprocessing. During the sweeps, data were acquired from analog 
transducers in the force balance and for tu nn el conditions. The model surface, test section 
sidewall, and wake rake pressures were acquired using an electronically scanned pressure 
system. The acquisition of these data was synchronized in time corresponding to the angle of 
attack sweep rate. Lift and pitching-moment coefficients were calculated from the force balance 
and from the integration of the surface pressure measurements. Good agreement between the 
integrated-pressure data and the force-balance data were obtained. In this report, the lift and 
pitching-moment data reported for the clean configuration were obtained from the surface 
pressures, while the data reported for the iced configuration were obtained from the force 
balance. Since the stall of the iced-airfoil configuration was characterized by unsteady flow, the 
available signal conditioning for the force-balance data allowed for more effective filtering of 
these unsteady effects. The drag coefficient was calculated from the wake pressures using 
standard momentum-deficit methods; these values are reported in this report for all 
configurations. The performance coefficients were corrected for wind tunnel wall effects using 
the methods of Allen and Vincenti [15]. The angle of attack sweeps were perfonned for a large 
range of Reynolds numbers and Mach numbers, as shown in table 1. The matrix was designed to 
isolate the independent effects of these parameters. Therefore, Reynolds number variations were 
performed at Mach numbers of 0.10 and 0.20, while a Mach number variation was performed at 
a nominal Reynolds number of 12.2xl0 6 . 
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Table 1. Matrix of Reynolds and Mach Number Conditions 


Reynolds 

Number 

Mach Number 

0.10 

0.20 

0.28 

4.7xl0 6 

X 



8.8xl0 6 

X 

X 


12.2xl0 6 

X 

X 

X 

16.0xl0 6 


X 



The experimental uncertainty in the performance coefficients was estimated using the methods of 
Kline and McClintock [16] and Coleman and Steele [17] for 20:1 odds. Table 2 lists these 
uncertainties for both integrated-pressure and force-balance measurements, before the wall 
corrections were applied. The values were calculated based upon the clean model configuration 
at Re = 8.1xl0 6 and M= 0.20. The absolute uncertainties in table 2 are inversely proportional to 
the dynamic pressure (except for a). This condition was selected because it corresponds to the 
average dynamic pressure over the range of conditions (table 1). Therefore, conditions having 
lower dynamic pressure would have slightly larger uncertainties, while conditions with higher 
dynamic pressure would have slightly lower uncertainties. All uncertainties were acceptable for 
the purposes of this investigation. The relative uncertainty in C m (both pressure and balance) 
seems large for this example because of the small reference value. For cases where the C,„ 
values were larger, e.g., in the iced-airfoil case, the absolute uncertainty would be similar, 
therefore resulting in a lower relative uncertainty. This is also the case for the uncertainty in 
drag coefficient. Finally, several repeat runs were performed for both clean and iced 
configurations and these run-to-run variations in the coefficients were much smaller than the 
uncertainties listed in table 2. 

Table 2. Estimated Experimental Uncertainties for Measurements in the 

ONERA FI Wind Tunnel 


Aerodynamic 

Quantity 

Reference 

Value 

(degrees) 

Absolute 

Uncertainty 

(degrees) 

Relative 

Uncertainty 

(%) 

a 

8.01 

±0.02 

±0.25 

Ci Balance 

1.095 

±0.010 

±0.93 

C m Balance 

-0.0144 

±0.00071 

±4.9 

Cp 

-1.057 

±0.032 

±3.05 

Ci Pressure 

1.096 

±0.0070 

±0.64 

C m Pressure 

-0.0148 

±0.0024 

±16.5 

C d Wake 

0.0086 

±0.00048 

±5.5 
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2.2 UNIVERSITY OF ILLINOIS WIND TUNNEL AND MODEL. 


The subscale aerodynamic testing was conducted in the University of Illinois subsonic, low- 
turbulence, open-return wind tunnel, which had a 33.63-inch (853-mm) high, 48-inch (1219-mm) 
wide, and 96-inch (2438-mm) long test section. The tests were performed on an aluminum 
NACA 23012 airfoil model with a chord length of 18 inches (457.2 mm) and a span of 
33.563 inches (852.5 mm). The model had a primary chordwise row of 73 pressure taps located 
at 51% span. Lift and pitching-moment coefficient data were acquired from a force balance and 
by integration of airfoil surface pressures measured by an electronically scanned pressure 
system. Excellent agreement between these methods was obtained for the clean model 
configuration. For the iced configurations, the ice simulations used in this study did not have 
pressure taps. This resulted in minor disagreement in the integrated coefficients with the force 
balance. Therefore, the lift and pitching-moment data shown in this report were obtained from 
the force balance. Using standard momentum-deficit methods, the drag coefficient was 
computed from total pressure measurements collected by a traversable wake rake. This wake 
rake is shown installed behind the airfoil model in figure 3. Busch [18] and Blumenthal [19] 
describe the experimental set-up in greater detail. Experimental uncertainties were calculated 
using the methods of Kline and McClintock [16] and Coleman and Steele [17], and a summary 
of these uncertainties is given in table 3. While the relative uncertainty of C m appears to be large 
because the reference value is small, the absolute uncertainty is reasonable. The angle of attack, 
lift, pitching-moment, and drag coefficients were corrected for wind tunnel wall boundary effects 
using the methods of Allen and Vincenti [15]. All data were collected at a Reynolds number of 
1.8xl0 6 and a Mach number of 0.18, unless noted otherwise. 



Figure 3. Photograph of Subscale NACA 23012 Airfoil Model Installed in University of Illinois 

Wind Tunnel Test Section 
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Table 3. Estimated Experimental Uncertainties for Measurements in the University of 

Illinois Wind Tunnel 


Aerodynamic 

Quantity 

Reference 

Value 

(degrees) 

Absolute 

Uncertainty 

(degrees) 

Relative 

Uncertainty 

(%) 

a 

4.13 

±0.02 

±0.48 

Ci Balance 

0.556 

±0.00086 

±0.16 

C m Balance 

-0.0015 

±0.00027 

±17.4 

Cp 

-0.962 

±0.0045 

±0.47 

C d Wake 

0.0071 

±0.00014 

±1.9 


2.3 FULL-SCALE MODEL ICE-SHAPE SIMULATION METHODS . 

The runback ice simulation tested on the full-scale model was a high-fidelity ice casting that was 
adapted from an earlier icing test conducted at the NASA IRT [20 and 21]. The runback ice 
accretion designated NG067 1 was considered to be representative of flight in holding conditions 
for a full-scale aircraft equipped with a bleed-air, thermal ice-protection system. A mold and 
casting of this ice accretion were produced and pictures of the completed casting are shown in 
figure 4. A tracing of the runback ice accretion NG0662 from a companion run is shown in 
figure 5. All test conditions between NG0671 and NG0662 were identical, so the tracing shown 
in figure 5 represents a cross section of the runback ice accretion shown in figure 4. Due to the 
operation of the thermal ice protection system, the leading edge is completely free of ice up to 
about x/c = 0.13 on the upper surface and x/c = 0.15 on the lower surface. Characteristics of the 
upper-surface ridge are that it was continuous along the span, had significant chordwise extent, 
and had some variation in height along the span. In contrast, the lower-surface ridge had a 
significant amount of spanwise variation with the largest ice elements being isolated, three- 
dimensional features, instead of fonning a continuous ridge. As discussed in section 1, the three- 
dimensional nature of this type of accretion requires a high-fidelity simulation to obtain accurate 
aerodynamic results. 

Since the IRT model used to generate ice accretion NG067 1 was not identical to the full-scale 
NACA 23012 model used for the aerodynamic tests, the ice casting was modified to the latter 
geometry. The castings of the ice ridges shown in figure 4 were extracted from the IRT model 
geometry and applied to a leading-edge section corresponding to the full-scale NACA 23012 
aerodynamic model. This representation of the runback ice accretion was then used to generate a 
series of ice-casting sections used for the aerodynamic tests through a process described in detail 
by Broeren, et al. [22], The completed installation is shown in figure 6. A tracing of the 
completed ice shape on the NACA 23012 airfoil leading edge along with the pressure tap 
locations is given in figure 7. This tracing does not show the height profile of the large, isolated 
ice elements on the lower surface due to a spanwise shift in the tracing location from that shown 
in figure 5. The tracing in figure 7 only extends to x/c = 0.20, since this was the physical limit of 
the ice-casting section applied to the full-scale NACA 23012 aerodynamic model. Any runback 
ice downstream of x/c = 0.20 was not present in the casting simulation used for the aerodynamic 
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tests. As shown in figure 5, this amounted to a very small amount of ice and was considered 
negligible in terms of its attendant aerodynamic effect. 



Figure 4. Photograph of NG0671 Runback Ice Accretion Casting 
(Left — Upper-Surface Ridge, Right — Lower-Surface Ridge) 



Figure 5. Tracing of NG0662 Runback Ice Accretion Formed at Conditions Identical to the 

NG0671 Casting Shown in Figure 4 
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Figure 6. Completed Installation of Runback Ridge Simulation NG067 1 on the Leading Edge of 

the Full-Scale NACA 23012 Model 



Figure 7. Tracing of Runback Ice Casting NG0671 on the NACA 23012 Leading Edge 
(Open circles indicate pressure orifice locations.) 

2.4 SUBSCALE MODEL ICE-SHAPE SIMULATION METHODS . 

The subscale NACA 23012 model was tested with several variations of NG0671 runback-ridge 
simulations. These simulations were constructed using commercial off-the-shelf materials 
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because the runback accretion was too small (upper surface) and too three-dimensional (lower 
surface) to model with rapid-prototyping methods such as stereolithography. A “build-up” 
approach was employed in which individual ice features were added to the upper and lower 
surface incrementally. This approach yielded an understanding of the aerodynamic effects of 
individual simulation features. The simulation methods and terminology follow those 
established in closely related research [1, 12, 18, 23, and 24], Table 4 summarizes the individual 
features that were combined to simulate the NG0671 runback-ridge ice accretion. The “SG” 
notation refers to “simple geometry” which, in this case, means that the runback ridge was 
simulated with a two-dimensional rectangular shape. This was attached to the model surface at 
x/c = 0.13, identical to the full-scale ice accretion location. The height of the upper-surface ridge 
was k/c = 0.0028, identical to the maximum height of the upper-surface tracing in figure 7. This 
was a geometric scaling of the traced cross-section. The actual ice accretion did have some 
variation in height along the span that was not modeled with the simple-geometry simulation. 
Following the “build-up” approach, carborundum roughness elements with height k/c = 0.0008 
were applied to the top of the simple-geometry ridge (SG+R-US in table 4) to model roughness 
and height variations in the actual accretion. These simulations are shown in figure 8. 
Analogous simple-geometry simulations were also developed and tested on the model lower 
surface. In this case, using a two-dimensional simulation was difficult to define due to the highly 
three-dimensional nature of the lower-surface ice accretion, thus it required engineering 
judgment. The height k/c = 0.0039 was selected for the lower-surface SG ridge based upon 
matching the combined frontal area of large isolated ice elements. The sensitivity of the lower- 
surface ice accretion simulation was detennined by developing a higher-fidelity, three- 
dimensional simulation. The individual isolated ice elements were geometrically scaled and 
attached to a thin substrate on the lower surface. In addition, simulated frozen rivulets were also 
modeled. These features are shown in figure 9. Aerodynamic tests were performed with various 
combinations of upper- and lower-surface simulations to quantify their effects. 

Table 4. Summary of Runback-Ridge Ice Simulations 


Simulation 

Description 

SG-US 

Rectangular ridge on upper surface, with height k/c = 0.0028 at x/c = 0.13, 
chordwise extent A x/c = 0.047 

SG-LS 

Rectangular ridge on lower surface, with height k/c = 0.0039 at x/c = 0.15, 
chordwise extent A x/c = 0.042 

SG+R-US 

Roughness elements, height k/c = 0.0008 applied to upper-surface ridge 

SG+R-LS 

Roughness elements, height k/c = 0.0008 applied to lower-surface ridge 

SG+Riv-US 

Simulated frozen rivulets, height k/c = 0.0017 and density of 18 rivulets per 
inch-span applied to upper-surface ridge 

3D-LS 

Three-dimensional simulation of lower-surface ridge ice elements, height 
varied between k/c = 0.0022 and 0.0106 

3D+Riv-LS 

Simulated frozen rivulets, height k/c = 0.0017 and density of 18 rivulets per 
inch-span applied to three-dimensional, lower-surface ice simulation 
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Figure 8. Photographs of Subscale Model Upper-Surface Ridge Simulations 
(Left — Simple-Geometry Ridge (SG-US), Right — Simple-Geometry Ridge With 

Applied Roughness (SG+R-US)) 




Figure 9. Photographs of Subscale Model Ridge Simulations (Left — Upper-Surface, Simple- 
Geometry Ridge With Simulated Frozen Rivulets (SG+Riv-US), Right — Lower-Surface, Three- 
Dimensional Simulation With Simulated Frozen Rivulets (3D+Riv-LS)) 
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3. RESULTS AND DISCUSSION. 


3.1 FULL-SCALE RESULTS . 

The effect of the high-fidelity NG067 1 runback ice simulation on the performance of the NACA 
23012 airfoil is illustrated in figure 10 for Re = 16.0 X 10 (1 and M = 0.20. Broeren, et al. [22] 
provide a detailed analysis of the clean NACA 23012 performance data establishing its validity 
compared to archival and computational data. For this condition, the maximum lift coefficient of 
the clean NACA 23012 was 1.82 at an angle of attack of 18.1°. The abrupt loss of lift at stall is 
indicative of leading-edge stall. For this stall type, boundary-layer separation occurs near the 
leading edge without subsequent reattachment, resulting in separated flow over the airfoil and the 
significant decrease in lift [25]. The pitching-moment coefficient was nearly independent of 
angle of attack up until about a = 10°. The effect of the NG0671 ice simulation was to reduce 
the lift-curve slope and maximum lift coefficient. The ice simulation produced angle of attack 
dependence of the pitching-moment coefficient at much lower angle of attack, commencing at 
approximately a = 4°. The presence of the runback ice simulation resulted in a greater than two- 
fold increase in drag coefficient over the angle of attack range shown in figure 10. The 
maximum lift coefficient was reduced to 1.51 at a sta ii = 15.0° in addition to a change in the 
character of the stall. Analysis of the surface pressure data and mini-tuft flow visualization 
indicated significant trailing-edge separation at maximum lift coefficient. However, the decrease 
in lift for a >16° appeared to be due to flow separation on the entire upper surface originating 
from the leading edge or from the spanwise ridge. Therefore, the iced-airfoil stall type was 
classified as a combination of trailing-edge and leading-edge stall. 

While these are significant perfonnance effects, they are not as severe as those found for other 
ice shapes tested on the full-scale NACA 23012 model. Plotted in figure 10 for comparison are 
data for three other ice-casting configurations taken from Broeren, et al. [22]. Tracings of these 
ice shapes are shown in figure 11. EG 1162 and EG 1164 were leading-edge ice shapes, while 
EG 1 159 was a tall spanwise ridge. The performance data show the increasing severity for each 
configuration in terms of reduced Ci jmax and a sta ii, increased pitching-moment slope, and 
increased drag. It is interesting to note that the streamwise ice shape EG1 162 resulted in lower 
drag than the runback ridge NG0671 up to a ~ 9.5°, as the airfoil with the EG1162 ice shape 
began to stall. As described in detail by Broeren, et al. [22], the stall of this configuration 
maintained the character of the abrupt, leading-edge stall type of the clean airfoil. For the 
EG 1164 and EG 1159 configurations, the aerodynamics were governed primarily by a large 
upper-surface separation bubble resulting from the ice accretion. The separation bubble grew 
larger with increasing angle of attack precipitating the stall at a much lower angle of attack than 
for the clean airfoil. This combined with the physical size of the ice shape also contributed to 
the large increase in drag from the clean airfoil configuration. 
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Figure 10. Comparison of Clean and Iced NACA 23012 Airfoil Performance for Various 
High-Fidelity, Ice-Casting Simulations at Re = 16.0 X 10 6 andM= 0.20 





Figure 11. Various Ice-casting simulations Tested by Broeren, et al. on the Full-Scale 
NACA 23012 Airfoil Model (Open circles indicate pressure orifice locations.) 
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The aerodynamic effect of these ice shapes is compared further in terms of surface pressure 
distribution in figure 12. Here, the data are split into two separate plots for clarity. First, 
consider the comparison of the NG0671 runback ridge, EG 1 162 streamwise ice, and clean airfoil 
configurations. At this angle of attack, there was very little difference in surface pressure 
between the clean airfoil and the EG 1162 configuration. The streamwise ice shape, being 
conformal to the NACA 23012 leading edge, does not result in large-scale boundary-layer 
separation or other significant redistribution of static pressure for this angle of attack. This 
explains, in part, the minimal increase in drag from the clean configuration relative to the other 
iced configurations for this angle of attack. The surface pressure distribution of the NG0671 
configuration clearly reveals the effect of the ice shape on both the upper and lower surface. At 
this positive angle of attack, the effect on the upper surface was much more pronounced. There 
was a large change in pressure, from C p = -0.4 to -2.7 at the forward face of the ridge. These 
measurements capture the flow deceleration immediately forward of the ridge followed by the 
flow acceleration over the top of the ridge. This was then followed by a rapid pressure recovery 
with the static pressure reconforming to the clean airfoil values of x/c >0.18. Surface-oil flow 
visualization was performed for this case by applying a thin film of oil treated with fluorescent 
dye to the surface of the model. The model angle of attack was set and the tunnel was operated 
allowing the oil to flow on the surface until a steady-state condition was achieved. These 
visualizations revealed the presence of a small separation bubble immediately aft of the ridge 
with subsequent reattachment within only a few percent chord. An important characteristic to 
note is that the pressure distribution was only altered from the clean configuration in the 
immediate vicinity of the runback ridge. Furthennore, the effect on the surface pressure was 
somewhat symmetric in that the effect on the integrated lift and pitching moment was small. The 
effect on pitching moment was also mitigated by the location of the ridge being close to the 
moment center at x/c = 0.25. Clearly the drag coefficient was increased by the forward face of 
the ridge and resulting separation. This explains, in part, why the drag coefficient was higher 
than for the EG1 162 configuration. There was also additional drag due to the ice elements on the 
lower surface of the NG067 1 runback-ridge configuration. 

Also plotted in figure 12 are surface pressure distributions for the EG1164 and EG1159 iced- 
airfoil configurations. For these cases, there was a large deviation from the clean airfoil owing to 
the large extent of separated flow. For the EG1 164 shape, the flow separated near the tip of the 
horn resulting in the region of nearly constant pressure from x/c = -0.02 to 0.03 on the upper 
surface. There was significant pressure recovery downstream of x/c = 0.03, but C p did not 
approach that of the clean airfoil until x/c = 0.20. For the EG1159 spanwise ridge shape, the 
separated flow region was much larger as indicated by the region of nearly constant pressure 
from x/c = 0.04 to 0.30 on the upper surface. This region was followed by a very gradual 
pressure recovery. Surface-oil flow visualization performed at this angle of attack indicated a 
time-averaged separation bubble reattachment zone from x/c = 0.64 to 0.68. The pressure data 
for this case also reveal the large effect of the lower-surface ridge. It also had a small separation 
bubble associated with it. The large extent of separated flow for the EG 1159 configuration is 
consistent with the large degradation in the performance coefficients in figure 10. 

The angle of attack evolution of surface pressure distribution is shown in detail in figure 13 for 
the NG0671 configuration. The iced-airfoil pressure distribution shown in figure 10 for a = 4.2° 
was maintained for all angles of attack prior to stall. That is, there was only a significant 
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deviation from the clean pressure distribution in the immediate vicinity of the upper and lower 
surface ridges. The comparison plot at a = 16.0° shows the deviation of the iced-airfoil C p 
toward the trailing-edge, indicating the presence of boundary-layer separation. Increasing the 
angle of attack to 17.4° resulted in complete separation of the upper-surface boundary layer, as 
reflected in figure 13. 




Figure 12. Comparison of Clean and Iced NACA 23012 Airfoil Pressure Distribution for 
Various High-Fidelity, Ice-Casting Simulations at a = 4.2°, Re = 16.0x1 0 6 , and M= 0.20 
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Figure 13. Comparison of Clean and Iced NACA 23012 Airfoil Pressure Distribution for 
Various Angles of Attack at Re = 16.0x1 0 6 and M= 0.20 

The pressurization capability of the ONERA FI wind tunnel was fully exploited to establish the 
effects of Reynolds and Mach number on the perfonnance of the NACA 23012 airfoil with the 
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NG0671 runback ridge. In general, changes in Reynolds number from 4.7 X 10 6 to 12.2 X 10 6 at 
M= 0.10 and from 8.9 X 10 6 to 1 6.0 X 1 0 6 at M = 0.20 had very little effect on the iced-airfoil 
performance, as illustrated in figures 14 and 15. Changes in Mach number from 0.10 to 0.28 at 
Re = 12.2 X 10 6 had slightly more effect on the iced-airfoil performance, with a small increase in 
drag coefficient over this range, as shown in figure 16. These results are consistent with 
previous aerodynamic studies of iced airfoils [8, 22, and 26-29]. 



e — 

- Re = 12.2xl0 6 , Af = 0. 10 

A 

- Re= 12.3xl0 6 , Vf = 0.20 

e — 

- Re = 12.2xl0 6 , iW = 0.28 



Figure 14. Effect of Reynolds Number on the Perfonnance of the NACA 23012 Airfoil With the 

NG0671 Ice Casting at M= 0.10 
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- Re = 12.3xl0 6 , M= 0.20 

B 

- Re = 16.0xl0 6 , M= 0.19 



Figure 15. Effect of Reynolds Number on the Perfonnance of the NACA 23012 Airfoil With the 

NG067 1 Ice Casting at M = 0.20 
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Figure 16. Effect of Mach Number on the Performance of the NACA 23012 Airfoil With the 

NG0671 Ice Casting at Re - 12.2><10 6 

The effect of Reynolds and Mach number on iced-airfoil performance can be further summarized 
through analysis of the maximum lift coefficient. This parameter is plotted for various 
configurations against Reynolds number in figure 17. These data, along with the data from 
Broeren, et al. [22], were all acquired in the ONERA FI wind tunnel. The data from Broeren, 
et al. [8] were acquired in the NASA Langley Low-Turbulence Pressure Tunnel. Since both 
facilities were pressurized wind tunnels, the Mach number was held constant at the values 
specified. The comparison in Ci, max for the clean NACA 23012 airfoil is considered to be good, 
given the potential differences in models, model scale, installation, facilities, as well as 
experimental uncertainty. The Ci jmax increased from 1.50 at Re = 2.0 X 10 6 , M = 0.12 to 1.88 at 
Re = 12.3 X 10 6 , M = 0.10. As shown in figure 17, there was much less Reynolds number 
dependence of Ci jmax for the iced-airfoil configurations. There was some variation in Ci, max with 
Reynolds number for the NG0671 runback ridge, but no clear trend. More data are needed to 
establish what Reynolds number dependence there may be, particularly for Re < 4.7 X 10 6 . 
Included for comparison to the NG0671 runback ridge are the four ice shapes from figures 10-12 
and a forward-facing, quarter-round spanwise ridge ice shape. The quarter-round shape had 
height k/c = 0.0139 and was tested at two chordwise locations on the airfoil ( x/c = 0.02 and 
0.10). The data for the leading-edge ice shapes EG1 162 and EG 1 164 show virtually no change 
in Ci max over the Re range tested. The data for the EG 1159 configuration indicate a small 
decrease in Ci, max for Re = 4.6 X 10 6 to 12.2 X 10 6 at M = 0.10. The quarter-round cases show even 
less variation with no significant change in Ci, max versus Re from 2.0 X 10 6 to 1 0.5 x 1 0 6 for these 
tall spanwise ridge shapes. The data in figure 17 show that Mach number is an equally, if not 
more, significant parameter for spanwise ridge ice shapes. 

The effect of Mach number on maximum lift coefficient is summarized in figure 1 8 for the same 
iced-airfoil configurations. The data show good agreement for the values of Ci, max on the clean 
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airfoil, particularly for M >0.20. Of the iced-airfoil configurations, the spanwise ridges had the 
most dependence on Mach number. For the NG0671 runback ridge, the Mach number 
dependence was similar to, if not more than, the clean airfoil. The Ci, max for the airfoil with the 
tall spanwise ridge shapes (EG 1159 and quarter round) also shows some Mach number 
dependence. Iced-airfoil Mach number effects have been addressed in other studies [2 and 30]. 

Aside from the Reynolds and Mach number effects, a more significant effect for the tall ridge 
shapes is the upper-surface location of the ridges. Close inspection of figure 1 1 yields a ridge 
height, k/c = 0.013 at x/c = 0.05 in chordwise location. This is very similar in height to the 
k/c = 0.0139 quarter-round shape data in figures 17 and 18. Lee and Bragg [7] researched the 
impact of ridge location and airfoil type on aerodynamic performance degradation. They found 
that Ci, max was a strong function of spanwise ridge location between x/c = 0.0 and 0.10 for this 
airfoil. This explains why there is such a large difference in Ci, max between the quarter-round 
shape at x/c = 0.02 and 0.10 and why the data for the EG 1 159 ridge fit neatly in between. 


Clean, M= 0.10, Ref. 22 

Clean, A/= 0.20, Ref. 22 

Clean, M = 0.12, Ref. 8 

NG0671 Runback Ridge, M= 0. 10 

NG0671 Runback Ridge, M= 0.20 

EG 1 162 Streamwise, M= 0.10, Ref. 22 

EG1 162 Streamwise, M= 0.20, Ref. 22 

EG 1 164 Horn, Af = 0.10, Ref. 22 

EG1 164 Horn, M= 0.20, Ref. 22 

EG1 159 Spanwise Ridge, M= 0.10, Ref. 22 

EG1 159 Spanwise Ridge, M= 0.20, Ref. 22 

Quarter Round at x/c = 0.02, M — 0. 12, Ref. 8 

Quarter Round at x/c = 0.02, M— 0.21, Ref. 8 

Quarter Round at x/c = 0. 10, M = 0. 12, Ref. 8 

Quarter Round at x/c = 0. 1 0, M = 0.2 1 , Ref. 8 
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Figure 17. Effect of Reynolds Number on Maximum Lift Coefficient for the NACA 23012 

Airfoil for Clean and Iced Configurations 
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Spanwise Ridge at x/c = 0. 10, Re = 7.5x 10 6 , Ref. 8 



Figure 18. Effect of Mach Number on Maximum Lift Coefficient for the NACA 23012 Airfoil 

for Clean and Iced Configurations 

3.2 SUBSCALE SIMULATION OF RUNBACK RIDGE AERODYNAMICS . 

Simulating the aerodynamics of a runback ridge on a quarter-scale model at low Reynolds 
number can be challenging because of Reynolds number effects for the clean airfoil. This 
conundrum is illustrated in figure 19, which shows that the iced-airfoil configuration at 
Re = 16.0x1 () 6 and M= 0.20 had a slightly higher lift-curve slope, Ci, max , and a sta u than the clean 
airfoil at Re = 1.8 x 10 6 and M = 0.18. The iced-airfoil configuration also had slightly less 
dependence of C m on angle of attack. Only the drag coefficient was significantly different from 
the low Reynolds number clean configuration. These data may explain the apparent anomalies 
in previous subscale runback-ridge aerodynamic studies. For example, Calay, et al. [4] observed 
small increases in Ci, max and a sta ii with a k/c = 0.0035 simple-geometry ridge on a NACA 0012 
airfoil at Re = 1.25xl0 6 and M = 0.08. Similarly, Papadakis and Gile-Laflin [5] also observed 
lift perfonnance increases for simple-geometry ridges on a NACA 63 A -213 airfoil at 
Re = 2.0xl0 6 and M = 0.17. Whalen, et al. [9 and 10] and Lee, et al. [13] are further examples. 
As shown in figure 19, the decrease in maximum lift coefficient and stalling angle of attack for 
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the clean airfoil from Re = 1 6.0 X 1 0 6 to Re = 1.8 X 10 6 was larger than the degradation due to the 
NG0671 ice-casting simulation at Re = 16.0x1 0 6 . Therefore, subscale simulations of the 
NG0671 shape at low Reynolds number must result in a slight increase in maximum lift and 
stalling angle relative to the clean airfoil to simulate the aerodynamic effect at higher Reynolds 
number. In the current study, the challenge was to identify the accuracy to which the 
aerodynamics of the NG0671 runback ridge could be simulated on the quarter-scale model at low 
Reynolds number. 



Figure 19. Comparison of C/, C m , and Cd vs a for the Clean Subscale Model and the 
Full-Scale NG067 1 Casting and Clean Configurations 


The present approach employed simple-geometry ridge simulations in various configurations. 
The initial perfonnance results are shown in figure 20 for three subscale configurations at 
f?e=1.8 x 10 6 and M = 0.18 compared to the full-scale NG0671 ice-casting and clean 
configurations at Re = 1 6.0 X 1 0 6 and M = 0.20. The data for the SG-US configuration show the 
effect of the k/c = 0.0028 upper-surface, simple-geometry ridge. The lift-curve slope for this 
configuration was slightly lower than the airfoil with the NG067 1 casting, but the Ci, max value of 
1.50 and a sta ii value of 15.4° compared very favorably with the full-scale data. The leading-edge 
stall characteristics of the clean airfoil persisted in the subscale iced-airfoil case in contrast to the 
full-scale casting simulation that had characteristics of trailing-edge stall. Despite the good 
comparison in Ci, max and a sta ii, the character of the iced-airfoil stall was not adequately simulated 
with the simple-geometry ridge on the upper surface. It is interesting to note, however, that the 
SG-US simulation slightly increased both Ci, max and a sta n relative to the clean airfoil values at the 
same Reynolds number and Mach number; thus acting like a low Reynolds number lift- 
enhancement device. Figure 20 also shows that the SG-US subscale simulation caused a 
stronger dependence of the pitching-moment coefficient on angle of attack than for the full-scale, 
iced airfoil. The effect on drag coefficient was larger at higher angle of attack, since the SG-US 
simulation was located on the upper surface. The addition of the lower-surface ridge (SG-US & 
SG-LS) improved the drag coefficient comparison while having very little effect on the stalling 
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characteristics. This superposition of iced-airfoil effects has been observed in other studies of 
horn-type ice shapes [2 and 31] and spanwise ridge shapes [32], Therefore, it was expected that 
the addition of the lower-surface ridge would primarily affect the drag, with little or no effect on 
lift and pitching moment at higher angle of attack. 



Figure 20. Comparison of C/, C m , and C'd vs a for Various Subscale Runback-Ridge Simulations 
and the Full-Scale NG067 1 Casting and Clean Configurations 

The effect of added roughness was determined via carborundum grains with k/c = 0.0008 applied 
to each upper- and lower-surface simple-geometry ridges. The lift data show that this had a very 
significant effect on the maximum lift coefficient and stalling angle of attack such that the 
comparison with the full-scale casting configuration was poor. These results are consistent with 
the conclusions of Calay, et al. [4] and Whalen, et al. [9], who noted that the stalling 
characteristics of the airfoil with simulated runback ridges can be very sensitive to the height or 
geometry of the ridge. As expected, the addition of the roughness on the lower-surface ridge 
increased the drag coefficient for this configuration at lower angle of attack, thus improving the 
drag comparison with the full-scale casting. However, the large deviations beginning at 
approximately 2° were due to the roughness applied to the upper-surface ridge and were 
consistent with the lower maximum lift coefficient and stalling angle. 

The differences in the aerodynamic effect of these simulations are further illustrated in the 
surface pressure distributions. Figure 21 compares the pressure coefficients for the upper- and 
lower-surface, simple-geometry ridges with and without applied roughness to the full-scale 
NG0671 casting and clean configurations at a matched angle of attack of 12.4°. The three 
runback-ridge simulations share some similarities. There were no significant differences 
between the clean and iced pressure distributions, except in the vicinity of the upper-surface 
ridge. For all iced configurations, the presence of the upper-surface ridge caused a small 
reduction in the magnitude of the suction peak and lower suction pressures downstream on the 
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upper-surface relative to the clean airfoil. The pressure coefficients near the trailing edge, 
particularly for the subscale ice simulations, were also lower, possibility indicating local 
boundary-layer separation at this angle of attack. The pressure distribution for the full-scale 
NG0671 casting configuration compared very favorably with the upper- and lower-surface, 
simple-geometry ridge simulation (SG-US & SG-LS). The large suction peak near the forward 
face of the ridge for the casting configuration was not measured for the subscale simulations 
since the subscale simulations were not instrumented with pressure taps. The favorable 
comparison between the pressure distributions for the NG067 1 casting and the subscale simple- 
geometry ridges (without added roughness) further confirms that there is proper simulation of the 
aerodynamics. The pressure distribution for the simple-geometry ridges with added roughness 
shows a deviation in C p just downstream of x/c = 0.17 on the upper surface. This deviation 
indicates the presence of a separation bubble downstream of the ridge with roughness. Since this 
separation bubble had a noticeable effect on the pressure distribution, it was likely larger than the 
separation bubble documented for the full-scale NG067 1 configuration. This larger separation 
bubble contributed to the stall at lower angle of attack compared to the no-roughness case, as 
well as the increased drag at higher angle of attack. 



Figure 2 1 . Comparison of Surface Pressure Distribution for Selected Subscale Runback-Ridge 
Simulations (at Re = 1.8 X 10 6 , M= 0.18) and the Full-Scale NG0671 Casting and Clean 
Configurations (at Re= 16.0 X 10 6 , M= 0.20), all at a = 12.4° 

Drag coefficient comparisons are more difficult to evaluate since there is usually not a single 
value of importance, as there is with maximum lift coefficient. A useful metric for comparing Cd 
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at multiple angles of attack is the percent root mean square (rms) difference in Q over an 
appropriate angle of attack range [12]. C/ varied linearly with a over the angle-of-attack range 
used in this study. The value A Cd,r ms is a percentage and is computed by determining the rms of 
the percent difference in Cd between the full-scale casting and subscale simulation at each angle 
of attack in the linear lift-coefficient range (a total of N angles of attack): 


AC 


d.rms 



( 1 ) 


While this parameter is useful for making comparisons over an angle of attack range, a 
disadvantage is that it does not give any indication as to whether the Cd of a given subscale 
simulation is higher or lower than that of the full-scale casting. For this study, the selected angle 
of attack range was -4° to 12°. Figure 20 indicates that the SG-US & SG-LS subscale simulation 
had the best overall comparison in Cd to the airfoil with the NG067 1 casting. The A Cd,rms for this 
case was 13.0%. The other two subscale simulations had A Cd, rm s values that were nearly twice as 
large. For comparison, the A Cd,,-ms for the full-scale NG0671 casting relative to the clean airfoil 
was 162%, indicating the large overall increase in drag due to the ice shape. 

Since the addition of roughness to the upper-surface simple-geometry ridge resulted in a large 
deviation in Cd relative to the NG0671 casting configuration, data were acquired with 
roughnesses added to the lower-surface ridge only. These data, labeled SG-US & SG+R-LS, are 
shown in figure 22. As expected, the stalling characteristics were similar to the configuration 
with only the upper-surface simple-geometry ridge (figure 20), while the agreement in drag 
coefficient was also improved. The A Cd.rms value was reduced to 10.2% for this configuration. 
Attempts at further improving the aerodynamic simulation were not successful. A three- 
dimensional, lower-surface simulation was fabricated that closely modeled the individual ice 
elements of the NG0671 shape. This was tested along with the simple-geometry ridge on the 
upper surface (SG-US & 3D-LS in figure 22). As expected, the influence of this change was 
primarily observed in the drag coefficient for a <4°. The corresponding values of Cd were 
significantly higher than for the airfoil with the NG067 1 casting. This led to a A Cd,rms value of 
14.3%. Finally, simulated frozen rivulets were applied to the upper-surface simple-geometry 
ridge and to the lower-surface three-dimensional simulation (SG+Riv-US & 3D+Riv-LS in 
figure 22). This simulation was tested because it represents perhaps the highest-fidelity subscale 
simulation possible using commercial off-the-shelf materials. Such a simulation may be 
considered “the most representative” without a priori knowledge of the full-scale, iced-airfoil 
aerodynamics. As shown in figure 22, the effect on performance was very similar to the addition 
of roughness to the simple-geometry ridges. The simulated frozen rivulets on the upper surface 
caused a reduction in maximum lift coefficient and stalling angle. The drag at angles of attack 
larger than 2° was significantly increased. These results are consistent with Whalen, et al. [9], 
who noted that three-dimensional runback ridge simulations tended to yield larger perfonnance 
penalties than their two-dimensional counterparts. 
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Figure 22. Comparison of C/, C m , and C'd vs a for Various Subscale Runback-Ridge Simulations 
and the Full-Scale NG067 1 Casting and Clean Configurations 

The subscale simulation effectiveness is summarized for each configuration in table 5 where 
Ci,max, otstaii, and A Cd,rms were computed relative to the full-scale NG067 1 casting configuration. 
According to these parameters, the best simulation was the simple-geometry ridge on the upper 
and lower surface with roughness applied to the lower-surface ridge only. This adequately 
simulated the full-scale iced-airfoil aerodynamics, with the exception of the stall type. This 
simulation on the subscale model stalled from the leading edge, in contrast to the full-scale iced 
airfoil that had a more gradual trading-edge stall. However, none of the subscale simulations 
tested were able to reproduce the stall characteristics of the full-scale iced airfoil. In addition, 
the pitching-moment coefficient also exhibited greater angle of attack dependence for all of the 
subscale simulation configurations relative to the full-scale iced airfoil. As previously discussed, 
the addition of roughness to the upper-surface simple-geometry ridge resulted in larger 
performance degradations. Similar results, though less severe, were obtained when simulated 
frozen rivulets were applied to the upper-surface simple-geometry ridge. Variations in the 
lower-surface ridge simulations primarily affected drag coefficient, predominately at lower 
angles of attack. 
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Table 5. Summary of Subscale Model Simulation Effectiveness for the NG067 1 Ice Casting 


Simulation 

Simulation 
Cl, max - Casting 
Cl, max 

Simulation 
C^stall " Casting (Xstall 
(degrees) 

Simulation vs 
Casting AQ„ 
(%) 

SG-US 

-0.01 

0.4 

25.6 

SG-US & SG-LS 

0.00 

0.4 

13.0 

SG+R-US & SG+R-LS 

-0.22 

-1.6 

23.8 

SG-US & SG+R-LS 

0.01 

0.4 

10.2 

SG-US & 3D-LS 

0.01 

0.4 

14.3 

SG+Riv-US & SG+Riv-LS 

-0.10 

-0.6 

16.7 


These data support the findings of previous researchers in regard to the aerodynamic effects of 
short, spanwise ridge ice shapes. At the subscale model Reynolds number, the simulated 
runback ice yielded slight increases in maximum lift coefficient and stalling angle compared to 
the clean configuration. The full-scale data from this study have shown that this apparent 
anomaly is due to the reduction in clean airfoil lift perfonnance with Reynolds number. It is 
important to note that an exact scale model of the NG067 1 casting was not tested on the subscale 
model nor was an exact scale model of the simple-geometry ridges tested on the full-scale model. 
This represents a Reynolds number gap between Re = 1.8 X 10 6 and 4.7 X 10 6 and a knowledge gap 
in the understanding of runback-ridge aerodynamics. The current research corroborates previous 
findings that the geometric details of the subscale simulation can have a significant effect on the 
iced-airfoil aerodynamics. Design and fabrication of three-dimensional subscale simulations or 
adding roughness are consistent with common practices. The present results indicate that such 
simulations tested at low Reynolds number will likely result in conservative performance 
degradations. This, too, is not well understood and may indicate that there are significant 
Reynolds number effects between Re = 1.8 X 10 6 and 4.7 X 10 6 . At this time, it is not clear why a 
simple, two-dimensional, geometrically scaled ridge simulation is required to best represent the 
full-scale, high Reynolds number, iced-airfoil aerodynamics. 

The boundary-layer based scaling methods investigated by Whalen, et al. [9, 10, 21, and 32] 
appear to produce subscale model perfonnance that is unrealistically degraded from the full-scale 
case. In that research, the boundary-layer scaled ice shapes were nearly a factor of two larger 
than the geometrically scaled ice shapes. According to the present data, the larger boundary- 
layer-scaled simulations would result in unrealistically large performance penalties compared to 
the full-scale, high Reynolds number case. Simulation methods for Reynolds numbers lower 
than 1.8 X 10 6 may require alternatives to geometric scaling, as was the case for Lee, et al. [13]. 
Figure 23 shows the change in perfonnance for the best simple-geometry simulation of the 
NG0671 casting with Reynolds and Mach number decreased to 1 .0 X 1 0 6 and 0.10, respectively. 
For the lower Reynolds number, the lift-curve slope approaching stall is reduced along with the 
maximum lift coefficient and stalling angle of attack. Interestingly, the pitching-moment and 
drag coefficients are largely unaffected by the decrease in Reynolds and Mach number. These 
data imply that an artificial ice shape smaller than the geometry-based scaling is required to 
adequately simulate the full-scale aerodynamics. The appropriate scaling length is unknown. 
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Figure 23. Effect of Reynolds and Mach Number on the Subscale Runback-Ridge Simulation 

Performance Effects 

3.3 AERODYNAMIC CLASSIFICATION OF SPANWISE RIDGE ICE . 

The data in figure 10 illustrate the contrasting aerodynamic effects between short spanwise 
ridges, as exemplified by the NG0671 shape, and tall spanwise ridges, as exemplified by the 
EG 1159 shape. Clearly, the aerodynamics are fundamentally different. The ice accretion 
classification developed by Bragg, et al. [2] for spanwise ridges was primarily based upon 
research related to tall ridges similar to the EG1159 shape. In this report, a subclassification is 
proposed that divides spanwise ridge ice into short and tall ridges. The definitions of these terms 
are not to be taken in regard to the physical height, but rather the associated effect on the flow 
field and pressure distribution. Specifically, the distinction between short and tall ridges lies in 
consideration of the separation bubble generated by the presence of the spanwise ridge and the 
resulting aerodynamics. 

The subclassification is analogous to the separation bubble nomenclature proposed by Tani [33] 
in 1964. In his review, Tani studied airfoil separation bubble characteristics. He defined the 
terms short and long, not based upon length, but upon the effect on pressure distribution. Tani 
wrote, “The presence of a long bubble makes the pressure distribution radically different from 
that in inviscid flow, with the result that the sharp suction peak near the leading edge is not 
realized. Instead, a suction plateau of a reduced level extends over the region occupied by the 
bubble length.” In contrast, a short bubble has a minimal effect on the pressure distribution and 
leading-edge suction peak. The distinction between short and tall spanwise ridge ice can be 
made in terms of the resulting separation bubble and its effect on the flow field as manifest in the 
pressure distribution. 
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Numerous investigations into the aerodynamics of tall spanwise ridge ice shapes have been 
conducted [6-8 and 34-36] and the results are very briefly summarized here. As previously 
mentioned, the salient feature of tall spanwise ridge flow fields is the separation bubble that 
forms downstream of the ice shape. The resulting redistribution of surface pressure is significant 
for the EG 1159 shape, as illustrated in figure 12. The iced-airfoil pressure distribution bears 
almost no similarity to the clean airfoil pressure distribution at the same angle of attack. Lee and 
Bragg [6] documented the growth of the upper-surface separation bubble with angle of attack for 
the tall, forward-facing, quarter-round, spanwise ridge shape using surface-oil flow visualization. 
The k/c = 0.0139 ice shape was located at x/c = 0.10 on the upper surface of the NACA 23012m 
airfoil and tested at Re - 1.8x1 0 6 and M= 0.18. The bubble reattachment region was centered at 
about x/c = 0.42 at a = 0°. The bubble grew rapidly with increasing angle of attack, until the 
reattachment region reached the trailing edge at about a = 4°, precipitating the stall with 
maximum lift coefficient of about 0.25. As discussed at length by Bragg, et al. [2], these large 
separation bubbles are known to have large-scale unsteady characteristics. Since the surface-oil 
flow visualization technique is essentially a time-averaged method, only the region of mean 
reattachment was indicated. Also investigated were the effect of ridge height, cross-sectional 
geometry, location, and airfoil geometry. The height of the spanwise ridge was found to have a 
strong influence on the iced-airfoil aerodynamics as well as the chordwise location of the ridge 
for certain airfoils. Cross-sectional geometry was also found to be an important factor. As 
discussed in connection with figure 17, there were negligible Reynolds number effects on 
maximum lift coefficient at constant Mach number down to Re = 2.0 X 10 6 . Lee and Bragg 
[6 and 7] showed continued Reynolds number independence down to Re = l.OxlO 6 andM= 0.10 
for the k/c = 0.0139 quarter-round shape located at x/c = 0. 10. It is clear that tall spanwise ridge 
shapes have a large effect on the airfoil pressure distribution in terms of deviation from the clean 
configuration. This is manifested through a large separation bubble that grows in size (in the 
time-averaged sense) with small increases in angle of attack. Linally, the effect of Reynolds 
number on the flow field and performance has been shown to be negligible for the iced airfoil 
down to Re = 1 .Ox 1 0 6 . 

Short spanwise ridges exemplified by the NG0671 shape have different characteristics. As 
discussed in connection with figure 12 (also shown in figure 13), the effect of the simulated 
runback ridge on the pressure distribution was only a localized disturbance in the immediate 
vicinity of the ridge (both upper and lower surfaces for the NG0671 configuration); elsewhere, 
the surface pressure closely matched that of the clean airfoil. Like Tani’s short and long bubbles, 
this distinction — effect on pressure distribution — determines what is a short or tall spanwise 
ridge shape. The question then remains to determine how obvious this distinction may be for 
shapes of intermediate heights. Whalen [32] perfonned an extensive parametric study of two- 
dimensional, square-cylinder, simple-geometry simulations of spanwise ridges on a NACA 3415 
airfoil at Re = 1.8 x 10 6 and M= 0.18. The effect on pressure distribution for increasing ice-shape 
height for two different angles of attack is shown in figure 24. This airfoil model was equipped 
with a plain flap with a hinge located at x/c = 0.75. The flap was not deflected for these 
experiments and the lower-surface gap was sealed. The data for a = 10° show the increasing 
aerodynamic effect of the square -ridge shapes for increasing height. The large discontinuity in 
C p is observed at the ice-shape location of x/c = 0. 16 on the upper surface. Clearly, the pressure 
distribution for the smallest shape ( k/c = 0.0035) had features analogous to the NG0671 runback 
ridge. In the former case, the presence of the simple-geometry ridge caused a slight decrease in 
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the suction pressure from the leading-edge peak downstream to the ice shape (at x/c = 0.16) 
relative to the clean configuration. Increasing the angle of attack to 13° shows much better 
alignment with the clean surface pressures, except in the immediate vicinity of the ridge. 



Figure 24. Surface Pressure Effect of Square-Cylinder Spanwise Ridge Height at x/c = 0. 16 on 
the Upper Surface of the NACA 3415 Airfoil at Re = 1.8 X 10 6 , M— 0.18 [32] 
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This behavior is contrasted with the k/c = 0.0069 square shape. Even at the lower angle of attack 
of 10°, the effect of this simulated ice shape on the surface pressure was profound. It exhibited 
the classic characteristics of tall spanwise ridges. The pressure distribution for the iced airfoil 
had significant deviation from the clean airfoil configuration. The pressure plateau from x/c = 
0.17 to 0.24 and the subsequent gradual pressure recovery are indicative of a large separation 
bubble. This was confirmed via surface-oil flow visualization that indicated a time-averaged 
reattachment location at x/c = 0.45. The deviation of surface pressure from the clean 
configuration downstream of x/c = 0.65 on the upper surface was indicative of boundary-layer 
separation that was also confirmed in the surface-oil flow visualization. The increase in angle of 
attack to 16° for this configuration indicates a stalled flow condition. Based upon these 
observations, the distinction between the k/c = 0.0035 and 0.0069 shapes is clear. 

For the intermediate -height, simple-geometry ridges, the distinction between what would be 
considered tall or short may not be as obvious. However, further comparison of the pressure 
distributions in figure 24 reveals some important characteristics. For a = 10°, there are only 
minor differences in the pressure distributions for the k/c = 0.0044 and 0.0052 configurations. 
The biggest difference is that the airfoil with the k/c = 0.0044 shape attained a lower pressure just 
downstream of the ridge (x/c = 0.018) (having nearly the same value as for the k/c = 0.0035 
shape (C p ~ -2.9)). The k/c = 0.0044 shape also resulted in a steeper pressure recovery while the 
k/c = 0.0052 configuration exhibited a short pressure plateau region with more gradual pressure 
recovery. The latter is indicative of a large separation bubble. Therefore, based on the proposed 
definition, the k/c = 0.0044 ridge would be considered a short spanwise ridge, while the 
k/c = 0.0052 would be considered a tall spanwise ridge. This is confirmed in the pressure data 
for a = 13°. Here, the pressure distribution for the k/c = 0.0044 configuration conformed to the 
clean airfoil C p , except in the immediate vicinity of the ridge, while the pressure distribution for 
the k/c = 0.0052 configuration more clearly indicates the presence of a large separation bubble. 
Therefore, for the case of a spanwise ridge ice shape located at x/c = 0.16 on the upper surface of 
the NACA 3415 airfoil at Re = 1.8 X 10 6 and M = 0.18, the boundary between what is a short 
spanwise ridge and tall spanwise ridge lies somewhere between k/c = 0.0044 and 0.0056. Of 
course, this distinction is somewhat arbitrary. In reality, there is clearly an overlap region in the 
proposed subclassification as there is in any continuum. 

Whalen, et al. [9, 10, 21, and 32] further investigated the flow fields for the k/c = 0.0035 and 
0.0069 spanwise ridges using surface-oil flow visualizations and hot-wire, boundary-layer 
measurements. Consistent with previous findings for tall spanwise ridges was the presence of a 
large separation bubble that grew rapidly with angle of attack for the k/c = 0.0069 ridge. The 
flow field for the airfoil with the k/c = 0.0035 spanwise ridge was fundamentally different. 
Although there was a small separation bubble, the reattachment location was nearly constant 
over a large angle of attack range. Figure 25 is a plot of the bubble reattachment location versus 
angle of attack for three different configurations. Included with the square-cylinder shapes on 
the NACA 3415 airfoil are data from Fee and Bragg [6] for the k/c = 0.0139 quarter-round shape 
on the NACA 23012m airfoil. The data clearly illustrate the difference in angle of attack 
dependence of the separation bubble reattachment location. The error bars illustrate the size of 
the reattachment region due to the aforementioned unsteady characteristics associated with these 
large separation bubbles. In contrast to the tall shapes, the error bars for the short, k/c = 0.0035 
shape on the NACA 3415 airfoil illustrate the more steady character of the bubble. For the tall 
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ridge configurations, the stall was driven by the increasing size of the separation bubble (thin 
airfoil stall). But in the case of the short ridge configuration, the principle stall mechanism was 
from trailing-edge separation moving forward with increasing angle of attack (trailing-edge 
stall). 


■© NACA 3415, k/c = 0.0035 at x/c = 0.16, Ref. 32 

-A NACA 3415, k/c = 0.0069 at x/c = 0.16, Ref. 32 

-a NACA 23012m, k/c = 0.0139 at x/c = 0.10, Ref 2 



Figure 25. Comparison of Separation Bubble Reattachment Location Based on Surface-Oil Flow 
Visualization Results for Tall and Short Spanwise Ridges, all Data at Re= 1.8 x 1 0 6 , M = 0.18 

The differences in the separation bubble details are also illustrated in the single-component, hot- 
wire mean and rms velocity profiles in figures 26 and 27. The profiles were measured on the 
NACA 3415 airfoil upper surface for the k/c = 0.0035 and 0.0069 ridge shapes at a = 8°, 
corresponding to the surface-oil flow visualization data in figure 25. Since a single-component 
hot wire is not capable of indicating flow direction, the data close to the wall (n/k = 0 in figures 
24 and 25) must be carefully interpreted. First, for the k/c = 0.0035 case (figure 26), the mean 
and rms velocity profiles at x/c = 0.163 were measured immediately downstream of the ridge. 
The remaining profiles, at progressive downstream locations are consistent with an attached, 
developing turbulent boundary layer. For the larger ridge, the rms velocity profiles indicate a 
significant departure from the airfoil surface of the shear layer downstream of the ridge. The 
magnitude of the peak rms velocities are much larger for the larger ridge, further indicating the 
fundamentally unsteady nature of the separation region formed downstream of this shape. This 
is very different from the smaller ridge case where the shear layer remains close to the surface as 
indicated by the rms velocity profiles. Several of the mean velocity profiles downstream of the 
larger ridge have an inflection point consistent with separated flow in this region. These 
measurements, together with the surface pressure and surface-oil flow visualization data, indicate 
a large separation zone for the tall ridge and minimal flow disturbance for the short ridge. 
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Figure 26. Mean and rms Velocity Profiles Downstream of k/c = 0.0035 Square-Cylinder Ridge 


Shape Located at x/c = 0. 16 on the NACA 3415 Airfoil With a = 8°, Re 

and M= 0.18 [32] 
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Figure 27. Mean and rms Velocity Profiles Downstream of k/c = 0.0069 Square-Cylinder Ridge 
Shape Located at x/c = 0.16 on the NACA 3415 Airfoil With a = 8°, Re = 1.8x1 0 6 , 

and M= 0.18 [32] 
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The chordwise location of upper-surface spanwise ridges is known to have a significant effect on 
the aerodynamics. With respect to the present discussion, the chordwise location can influence 
what is classified as a short or tall spanwise ridge. Whalen [32] studied this effect for the 
k/c = 0.0035 square-cylinder shape on the NACA 3415 airfoil. The effect on pressure 
distribution is shown in figure 28 for two locations. The baseline location was x/c = 0.16, where 
the shape is clearly classified as “short.” Locating the same artificial ice shape at x/c = 0.06 had 
the effect of significantly altering the pressure distribution. Thus, the same geometry had 
characteristics of a tall spanwise ridge when located farther forward on the airfoil. The 
corresponding lift data indicated a significant reduction in maximum lift coefficient and stalling 
angle of attack relative to the clean configuration at the same Reynolds and Mach number. This 
stands in contrast to the lift enhancement due to the same artificial ice shape located at the 
downstream location, x/c = 0.16. This effect of chordwise location is consistent with the short 
spanwise ridge studies of Calay, et al. [4] and Papadakis, et al. [5], where lift enhancements were 
observed for downstream locations while the upstream locations resulted in significant 
degradations. The effect of chordwise location is driven by a number of factors, such as 
boundary-layer thic kn ess. The airfoil boundary layer is generally thinner near the leading edge 
(depending on the transition location), so that even relatively short ridges are many times larger 
than the local boundary-layer thickness. This and other factors can result in a larger separation 
bubble downstream of the ridge, which in turn results in a larger effect on the surface pressure 
distribution and larger perfonnance degradations. 



Figure 28. Surface Pressure Effect of k/c = 0.0035 Square-Cylinder Spanwise Ridge Chordwise 
Location on the Upper Surface of the NACA 3415 Airfoil at a = 12°, Re = 1.8><10 6 , 

andM= 0.18 [32] 
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4. CONCLUSIONS. 


This report presents the results of recent investigations into the aerodynamics of simulated 
runback ice accretion on airfoils. Aerodynamic testing was performed on a full-scale model 
using a high-fidelity, ice-casting simulation at near-flight Reynolds numbers (Re). In addition, 
follow-on subscale tests were conducted with low-fidelity simulations on a quarter-scale model 
at low Reynolds numbers. The high-fidelity, ice-casting simulation was fabricated from an ice 
accretion mold acquired during an icing tunnel test on a full-scale model with an anti-icing 
system. It was attached to the leading edge of a 72-inch (1828.8-mm) chord National Advisory 
Committee for Aeronautics (NACA) 23012 airfoil model. Aerodynamic perfonnance tests were 
conducted at the Office National d’Etudes et de Recherches Aerospatiales (French Aeronautics 
and Space Research Center) (ONERA) FI pressurized wind tunnel over a Re range of 4.7xl0 6 to 
16.0xl0 6 and a Mach number (M) range of 0.10 to 0.28. For Re = 16.0xl0 6 and M= 0.20, the 
simulated runback ice accretion on the airfoil decreased the maximum lift coefficient from 1.82 
to 1.51 and decreased the stalling angle of attack from 18.1° to 15.0°. The pitching-moment 
slope was also increased and the drag coefficient was increased by more than a factor of two. In 
general, the perfonnance effects were insensitive to Reynolds and Mach number changes over 
the range tested. While these aerodynamic penalties are significant, they are generally less than 
the penalties associated with full-scale, high-fidelity, leading-edge ice simulations, including 
roughness, on the NACA 23012 airfoil. 

The results from the full-scale tests were used to evaluate simulation methods for runback ice 
accretion. Aerodynamic tests were conducted on a quarter-scale NACA 23012 model (18-inch 
(457.2-mm) chord) at Re = 1.8xl0 6 and M = 0.18, using low-fidelity, geometrically scaled 
simulations of the full-scale casting. It was found that simple, two-dimensional simulations of 
the upper- and lower-surface runback ridges provided the best representation of the full-scale, 
high Reynolds number aerodynamics with the ice-casting simulation. This adequately simulated 
the full-scale, iced-airfoil aerodynamics, with the exception of the stall characteristics. The 
simple, two-dimensional simulation on the subscale model resulted in stalling from the leading 
edge, in contrast to the full-scale iced airfoil that had a more gradual trading-edge stall. 
However, none of the subscale simulations tested were able to reproduce the stall characteristics 
of the full-scale iced airfoil. The pitching-moment coefficient also exhibited greater angle of 
attack dependence for all of the subscale simulation configurations relative to the full-scale iced 
airfoil. Addition of grit roughness to the lower-surface simple-geometry ridge further improved 
the simulation accuracy in the drag coefficient. However, the addition of roughness to the upper- 
surface simple-geometry ridge resulted in larger performance degradations (in lift and drag) 
compared to the full-scale, high Reynolds number configuration. A higher-fidelity simulation of 
the runback ice accretion that included geometrically scaled, three-dimensional features also 
resulted in larger performance degradations than what was measured on the full-scale model. 

At this time, it is not clear why a simple, two-dimensional, geometrically scaled ridge simulation 
is required to best represent the full-scale, high Reynolds number iced-airfoil aerodynamics. The 
fact that the aerodynamic perfonnance of the subscale simulations is very sensitive to the 
geometric details, such as the addition of surface roughness, implies that Reynolds number 
effects may not be negligible for this type of ice accretion. An exact scale model of the NG067 1 
casting was not tested on the subscale model nor was an exact scale model of the simple- 
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geometry ridges tested on the full-scale model. This represents a Reynolds number gap between 
Re = 1.8 X 1 0 6 and 4.7 X 10 6 , and a knowledge gap in the understanding of runback-ridge 
aerodynamics. Therefore, more research specifically targeting this Reynolds number range is 
recommended. 

From both past and current research, the appropriate scaling length for the runback ice accretion 
is unclear. The previously proposed boundary-layer thickness length scale would result in 
subscale simulations that were too large and would thus result in larger perfonnance 
degradations compared to the full-scale, high Reynolds number configuration. Ideally, a length 
scale could be detennined that was not explicitly dependent upon height and location of the ice 
accretion. Such a length scale may be dependent on Reynolds number. More research is 
required to detennine the appropriate length scale for subscale ridge simulations, particularly for 
Reynolds numbers less than 1.8x1 0 6 . 

The results of this investigation were used develop a new subclassification of spanwise ridge ice 
that distinguishes between short and tall ridges. This subclassification is based upon the flow 
field and resulting aerodynamic characteristics, regardless of the physical size of the ridge and 
the ice accretion mechanism. Tall spanwise ridges have a profound effect on the airfoil flow 
field with a large (and often unsteady) separation bubble that grows rapidly with angle of attack, 
precipitating the stall at low-lift coefficient and angle of attack. Previous research has shown 
that the aerodynamic effect of tall spanwise ridge ice is insensitive to Reynolds number over a 
large range from Re > 1.0 X 10 6 . In contrast, short spanwise ridges are characterized by a small, 
stable separation bubble formed in the immediate vicinity of the ridge. This small separation 
zone results in a limited effect on the airfoil pressure distribution relative to the clean 
configuration. Significant differences in surface pressure are only observed at the ridge location. 
The separation bubble does not significantly increase in size with angle of attack. Thus, the stall 
of the iced airfoil is generally governed by trailing-edge separation moving forward with 
increasing angle of attack (trailing-edge stall). In some cases, the airfoil may stall from the 
leading edge or from the ridge itself. More research is needed to determine the appropriate 
simulation methods for short ridges, particularly for Reynolds numbers less than 1 . 8 X 10 6 . 
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